The Evolution of Cas A at Low Radio Frequencies 
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We have used archival 74 MHz VLA data spanning the last 15 years in combination 
with new data from the Long Wavelength Demonstrator Array (LWDA) and data from 
the literature covering the last 50 years to explore the evolution of Cas A at low radio 
frequencies. We find that the secular decrease of the flux density of Cas A at ~ 80 MHz 
is rather stable over five decades of time, decreasing at a rate of 0.7-0.8% yr~^. This is 
entirely consistent with previous estimates at frequencies as low as 38 MHz, indicating 
that the secular decrease is roughly the same at low frequencies, at least between 38 
and 80 MHz. We also find strong evidence for as many as four modes of flux density 
oscillation about the slower secular decrease with periods of 3.10 ± 0.02 yr, 5.1 ± 0.3 yr, 
9.0 lb 0.2 yr, and 24 it 2 yr. These are also consistent with fluctuations seen previously 
to occur on scales of a few years. These results provide compelling motivation for a 
thorough low frequency monitoring campaign of Cas A to constrain the nature and 
physical origins of these fluctuations, and to be able to better predict the flux density of 
any given epoch so that it may be used as a reliable low frequency calibrator. 

Subject headings: catalogs — surveys — radio continuum: general 



1. Introduction 

Cassiopeia A, or "Cas A" is one of the brightest radio sources in the sky, making it one of the 
most well studied supernova remnants (SNRs). For some time, it has been known that the radio 
luminosity of Cas A is fadirig and that the rate at which it fades varies with frequency (see, e.g.. 



Hoebom k Shakeshaft l[l96ll:IScott et al. lll969l:[Baars et al. I ll977l: iRees II 19901 : lO'Sullivan Green 



199S 



Martirossvan et al. 



2002 



Vinvaikin II2006I ) . The canonical expres sion for the se c ular d ecrease 



of the flux density of Cas A as a func tion of frequency is that given by iBaars et al. I (119771 ). How- 



ever, recent work has shown that the iBaars et al. I (119771 ) equation may exaggerate the frequency 



dependence ( O'Sullivan Sz Green 19991 ) and at low (< 100 MHz) frequencies, the equation signifl- 



cantly overestimates the magnitude of the secular decrease (IRees Ill990l : iMartirossyan et al. II2002I : 



^Naval Research Laboratory, Code 7213, 4555 Overlook Avenue SW, Washington, DC 20375-5351 
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Vinvaikin 1120061 ) . In addition, it has also been demonstrate d that this decrease is not smooth, but 
i n fact contains significant variations (see, e.g., iRees Ill99d ). usuahy on time scales of a few years 
(jVinvaikin &: Razin II2005I ). Placing better and more current constraints on the magnitude and 
time dependence of both the secular decrease and smaller scale variations of the flux density of Cas 
A at low frequencies will aid in the establishment of the exact physical causes of these temporal 
and spectral variations. 



The implementation of the 74 MHz system for the NRAO Very Large Array (VLA: lKassim et al. 



20071 ) in the early 1990s provided a new means to monitor and explore the evolution of Cas A at 
low frequencies. New northern hemisphere instruments currently be ing developed such as the Low 
Frequency Array (LOFAR) and the Long Wavelength Array (LWA; iKassim et al. II2006I ) will allow 
for continued and better monitoring in the future. We have obtained archival 74 MHz VLA ob- 
servations as well as new 74 MHz observations with prototype LWA instrumentation of both Cas 
A and Cygnus A (or, "Cyg A") to provide a temporal baseline of more than 15 years with typical 
spacings of about 1-3 years. We require Cyg A observations bec ause it is a compar ably bright 



low frequency source with a flux density of 17.35 kJy at 74 MHz (iBaars et al. 1119771 ). Cyg A is 



also known to be non-variable because it is a radio galaxy whose flux density at low frequencies 
is dominated by extended emission spanning several kiloparsecs making any variability on human 
time-scales negligible. This makes Cyg A an ideal calibrator /reference source for our efforts. When 
combined with other results within the literature, these data have helped us achieve the goal of 
better constraining the evolution of Cas A at low frequencies up to the current epoch. In §2, we 
describe the new data and the measurements made with it. The characteristics of the temporal 
variations of the flux density of Cas A gleaned using these data are described in §3, and the impact 
this has on our understanding of the SNR and the low frequency flux density scale is discussed in 
^4. 



2. The Data 

All new, unpublished VLA data were selected from the NRAO Data Archive Syste Within 
the archive, observing runs that used the VLA and the 74 MHz or "4-band" system (exact central 
V = 73.8 MHz) which included observations of both Cas A and Cyg A were searched for. We 
identified seven different data sets from six observing programs with a sufficient number of scans 
of each source t o obtain a cahbrate d total flux density of Cas A using a "clean component" model 



of Cyg A on the lBaars et al. I (119771 ) scale (currently available at http://lwa.nrl.navy.mil/tutorial/) 
which approximates the 74 MHz intensity of Cyg A with a series of 697 delta functions (or, clean 
components) in the image plane. To increase the time and u-v coverage, the two data sets from the 
same observing program in 1991 (VLA project AK272) were combined into a single data set within 
AIPS. We also added to the archival data calibrated visibility data from the VLA Low Frequency 



^ ht t ps : / / archive . nrao . edu/ archive/ 
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Sky Survey (VLSS; VLA project AG735; ICohen et al. Il2007i . kindly provided by A. Cohen), of a 



field that contains Cas A. Unlike the other observations, Cas A was not in the center of the field 
of view for this VLSS field, and the flux density obtained from these data was corrected for the 
response of the primary beam of the VLA antennas. 

To obtain an integrated flux density of Cas A from each data set, the data were first cali- 
brated with standard AIPS routines using the aforementioned clean component model of Cyg A. 
Specifically, the AIPS routine CALIB was used with the Cyg A model to solve for the complex gain 
for each VLA antenna and each scan of Cyg A within time intervals of 1-3 minutes, depending on 
the data set. The CALIB task also flags time steps at which no good solution can be obtained. 
Following this, the gains determined by CALIB were applied to the scans of Cas A using the task 
CLCAL to interpolate onto the times of the Cas A scans. 

Following this, the visibility data were tapered with a Gaussian with a full width at half 
maximum (FWHM) in the u-v plane that corresponds to a Gaussian with a FWHM of 600 arcsec 
in the image plane (0.30 kilo- wavelengths) , roughly three times the angular diameter of Cas A. This 
tapering essentially takes Cas A from a relatively complex extended source and smoothes it to a 
single Gaussian source. This implies that any phase errors present in the Cyg A-based calibrations 
will have a negligible effect on the mean tapered amplitude as a function of u-v radius. 

Since Cas A is substantially brighter than any other source within the 74 MHz VLA primary 
beam (the typical total sky flux density at 74 MHz within the VLA primary beam is about 300 
Jy; A. Cohen, private communication, 2009) or any sources of RFI, no attempt was made to 
account for other sources within the field of view or to correct for the effects of RFI. Instead, the 
mean amplitude of the tapered data was computed within bins of u-v radius, taking the standard 
deviation within each bin to be the la error which includes the effects of both thermal noise. The 
radial amplitude plots, including the binned data, can be seen in Fig. [TJ A non-linear least-squares 
fitting routine was then used to fit a Gaussian function plus a constant to the binned data (see 
Fig. [T]for plots of these fits with the data), where the constant allows for the effects of noise, RFI, 
other sources in the field of view, etc., and the center of the Gassian was fixed at a u-v radius of 
zero. The peak amplitude of the fitted Gaussian (i.e., the amplitude at a u-v radius of zero) was 
then taken to be the total 74 MHz flux density of Cas A at that epoch, and the covariance matrix 
computed by the fitting routine for the free parameters was used to extract the la error in the flux 
density. 

We note, however, that this uncertainty may not accurately reflect the true uncertainty in the 
measurement of the absolute flux density since it is not clear whether or not the variance with 
each u-v radius bin contains a contribution from the uncertainty in the Cyg A-based calibration. 
Given the relative sophistication of the AIPS routines used, simply estimating the uncertainty in 
the measured Cyg A brightness and propagating this uncertainty through to our final flux density 
computation is not practical. However, since each Cas A scan is essentially calibrated separately 
[i.e., the calibration of a particular Cas A scan is based on the Cyg A scan(s) that is(are) nearest 
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in time], we can obtain a separate estimate of the Cas A flux density uncertainty via jackknife 
resampling that contains the effects of calibration uncertainty. We have done this for each data set 
by re-determining the flux density of Cas A using the same method as described above, each time 
excluding a different Cas A scan. The standard deviation among these re-determined flux densities 
subsequently contains a contribution from the calibration uncertainty since this uncertainty will 
manifest as an artificial variation in flux density from one scan to the next. For three of the data 
sets (AK272, AK332, and AG 735), these uncertainties were larger than those obtained from the 
Gaussian fits displayed in Fig. [T] and errors computed via jackknife resampling were subsequently 
adopted as the la uncertainties in the total flux densities of Cas A for these data sets. For all data 
sets, the flux densities and their final la uncertainties are listed in Table [1] and are plotted versus 
epoch in Fig. [2j 

To supplement these new data, w e have also compiled a set of low frequency data from the liter- 



ature (see below ). Data at 81.5 MHz (iHogbom &: Shakesha: 



1961 



Rees 



1991 



Scott et al. 
d: Ivinyaikin 



1968: 



2001 



Baars et al. 
~^) were ob- 



19771 ). 80 MHz (|Martirossvan et al. II2002I ). and 38 MHz 
tained and verified to be reliable measurements with no appreciable confusion issues. To add one 
additional epoch, we have obtained an estimate of the current flux density of Cas A from efforts 
to measure the beam pattern of antennas being developed for the LWA from September 2008 to 
Februar y 2009. Using two L WA prototype antenna s on the Long Wavelength Demonstrator Array 



(LWDA; I York et al. Il2007l i^. | Hartman et al. I (j2009l ) estimated that the flux density ratio of Cas A 
to Cyg A at 74 MHz over this time period was 1.04 it 0.05. We note that all of the se additional 
sources of data either give Cas A flux densities that are on the iBaars et al. I (119771 ) flux density 
scale, or they directly report Cas A to Cyg A flux density ratios and associated errors. Therefore, 
they can be used with our new VLA data to provide a consistent measure of the evolution of the 
flux density of Cas A. 



3. Temporal Variations 
3.1. Secular Evolution 



It has long been established tha t the flux density of Ca.s A has been and is decreasing with 



time at all radio frequencies (see, e.g., 'Hogbo m Sz Shakeshaft 



1961 



Scott et al. 



1977 



Walczowski fc Smith 



1985 



19691 : 



Baars et al. 



Rees nl99Qi ). This is confirmed by the data plotted in Fig.[2]where 



the 74 MHz flux density deceases by a factor of about 1.4 b etween 1991 a nd 20 06. Comparing 
the VLA 74 MHz flux densities to the flux densities used by iBaars et al. I (119771 ) to derive their 
frequency dependent expression for the secular decrease of Cas A' s flux density revea ls that the 
recent VLA- measured flux densities are r nuch h igher than what the lBaars et al. I (|l977l ) expression 



predicts. Specifically, the iBaars et al. I (jl977l ) expression estimates a 1.3% yr ^ decrease at 74 



^see the LWA memo series currently located at Iittp://www.ece.vt.edu/swe/lwa/| 
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MHz. In con t rast, when we combi ned the Cas A to Cyg A flux density ratios from the data used by 



Baars et al. I (| 19771 ) at 81.5 MHz (jHogbom fc Shakeshaft Ill96ll : IScott et al. Ill969l ') with our newly 



presented data, t he 74 MHz LWDA data, a nd the 80 MHz Byurakan Astrophysical Observatory 
(BAO) data from iMartirossvan et al. I (|2002l '). we found something significantly different. First, we 
scaled the Cas A to Cyga A ratios from the VLA and LWDA 74 MHz data and the 8 MHz BAO 
data t o the expected values at 81.5 MHz using the Cas A and Cyg A spectra from [Baars et alTI 
(jl977l ). We then performed a non- linear le ast-squares fi t of a secular decrease model to these 
scaled data and the 81.5 MHz data used by iBaars et al. I (119771 ) and found a smaller decrease of 
0.77 ± 0.02% yr-i (see the upper panel of Fig. E]). 



This is similar to the 0.8% yr~^ decrease at 38 MHz found using data compiled bv lRees I ( 1990l ). 
In fact, by simply scaling the VLA and LWDA 74 MHz data and the BAO 8 MH z data according 
to the shapes of the spectra for Cas A and Cyg A given bv iBaars et al. I (|l977l ) , we have found 
that both data sets are entirely consistent with the results presented by iRees I (j 19901 ). with the 
possi ble exception of t he VLA data from 1993. The same is true for additional 38 MHz available 
from Ivinyaikin I ( 2006 ). This is illustrate d in th e lower pan el of Fig. [3l where we have plotted the 
scaled 74 and 80 M Hz d ata as well a s the iRees I (|l990l ) and IVinvaikin I (|2006l ) data with both the 



Baars et al. I (jl977l ) and iRees I (|l990l ) predictions for the secular decrease of Cas A at 38 MHz. A 



new fit to all these data, whi ch is al s o plot ted in Fig. [3l gives a secular decrease of 0.84ifc0.05% yr ^, 
nearly identi c al to that from iRees I (jl990l ). It therefore appears that, as noted bv iRees I (jl990l ). the 



Baars et al. I (jl977l ) expression for the secular decrease of the brightness of Cas A at low frequencies 



does not hold for recent epochs. Additionally, the magnitude of this secular evolution appears to 
be independent of frequency, at least for 38 < < 80 MHz, and to not have significantly changed 
in the last ~ 50 years. 



3.2. Smaller Scale Variations 

While the (38 MHz) data plotted in the lower panel of Fig. [3] demonstrate a clear, consistent 
picture of the secular decrease of the flux density of Cas A, there appears to be some evidence 
of variations on shorter time scales. This is indicated by the scatter of the data about the 0.84% 
yr~^ decrease model which ranges from being within the la uncertainties to as much as 2-3cj. 
Specifically, the reduced (x^ P^^ degree of freedom) is 1.35. For uncorrelated, random Gaussian 
errors, this implies there is a roughly 50% chance that the simple secular decrease model explains 
the data by itself. Assuming these uncertainties are believable, this points to real deviations from 
the smooth secular decrease of Cas A's brightness w hich appear to occur on ti me scales on the 



order of a few years similar to the variations noted by IVinyaikin Sz Razin I (120051 ) . To explore this 
further, we have computed the discrete Fourier transform (DFT) of the 38 MHz data minus the 
0.84% yr~^ decrease model. The mean separation in time between the data points is 1.5 yr and we 
have therefore only computed the DFT for |z^| < 0.34 yr~^, the approximate Nyquist limit. The 
amplitude of the resulting DFT is plotted in the upper panel of Fig. U] and has 3-6 relatively large 
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peaks. To explore whether or not the smaller scale fluctuations could be approximated with a sum 
of regular oscillation modes, we used the properties of the DFT near the six largest peaks as initial 
guesses for a nonlinear least-squares fit to the data of sums of cosines. For models with sums of 
1, 2, 3, 4, 5, and 6 cosines, the reduced values (taking into account the two fewer degrees of 
freedom due to the two parameters of the fitted secular decrease model) were, 1.05, 0.777, 0.762, 
0.670, 0.687, and 0.787, respectively. 

From the x^-based analysis above, it appears that the smaller scale variations about the secular 
decrease model are well approximated by a sum of four cosines and that adding more cosines does 
not improve the fit. However, we must note that the quantity is only an adequate measure of the 
goodness of a fit if the errors used are reliable measurements of the uncertainties in the data. It is 
possible that the true uncertainties in several of the data points have been underestimated and that 
the spectrum seen in Fig. [4] is the result of random variations about the secular decrease of Cas A's 
flux density that happen to be reproducible with a sum of four cosines. To check this, we generated 
10^ simulated data sets, each time randomly drawing 37 numbers [i.e., the number of epochs for the 
comb inatio n of th e 74 MHz VLA and LWDA data, the 80 MHz BAO data, and the 38 MHz data 
from iRees I (jl990l ) and from IVinvaikin &: Razin I (|2005l )] from a normal distribution with a mean 
of zero and standard deviation equal to the 0.11 rms deviation of the data about the 0.84% yr~^ 
decrease model. An examination of the cumulative distribution of these deviations revealed that 
a normal distribution with a mean of zero and standard deviation of 0.11 approximates the data 
well. For each of the simulated data sets, the DFT was computed using the epochs of the real data. 
We then used the cumulative distribution of the amplitudes from all 10^ DFTs at each frequency 
to compute the 75%, 90%, and 95% confidence levels as functions of frequency. We have plotted 
these functions with the DFT of the data in Fig. HI From this plot, it can be seen that in fact, 
the four largest peaks have amplitudes above the 95% confidence level which we have highlighted 
with grey points. This confirms what the x^-based analysis implied, that a four-cosine model is a 
good representation of the real, smaller scale temporal variations in the brightness of Cas A. The 
best-fitting four-cosine model parameters and their estimated la errors are listed in Table [2j 

We have plotted the (38 MHz) data (minus the 0.84% yr~^ decrease model) in the lower panel 
of Fig. [4] along with the best fitting four-cosine model. We have also shaded the area corresponding 
to the la uncertainty in the model which we computed using the covariance matrix of the free 
parameters and the partial derivatives of the model function with respect to the free parameters 
at each epoch. The model appears to represent the data well with the median absolute deviation 
from the model fit being about O.Scr. It can also be seen from Fig. U] that the model approaches 
the largest deviant points at epochs of 1961.5, 1967, and 1993.5. We have also plotted the DFT 
of the model fit evaluated at the epochs of the data in the upper panel of Fig. [H The main peaks 
are reproduced well in both location and amplitude while only the locations of some of the smaller 
peaks are repoduced. The fact that the amplitudes of the smaller peaks for the DFT of the data 
are larger is likely due to the influence of noise as well as that of higher frequency oscillations which 
cannot be fully explored here given the typical sampling rate of our data, but which likely have an 
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effect at lower frequencies due to secondary "sine" peaks caused by the Fourier transform of the 54 
yr rectangular envelope of the data (i.e., the data spans epochs from 1955 to 2009). This points to 
a distinct possibility that Cas A has (at least) four regular modes of flux density oscillations at low 
radio frequencies with periods of 3.10 it 0.02 yr, 5.1 it 0.3 yr, 9.0 it 0.2 yr, and 24 it 2 according to 
the best fitting four-cosine model. 



4. Discussion 

By combining archival 74 MHz VLA data and the new LWDA data with data from the litera- 
ture, we have successf ully confirmed that the secular decrease of Cas A at low frequencies is in fact 



not as extreme as the iBaars et al. I (jl977l ) results imply. We have also shown with this data that 



the low frequency flux density of Cas A oscillates on tim e scales of ^5-10 ye ars with possibly four 



or more modes. These were apparently the cause of the iBaars et al. I (|l977l ) misestimation of the 
magnitude of the secular decrease of Cas A as they did not have data that spanned a long enough 
time range and unfortunately fit a model to data that included a periodic "dip" in the flux density 
(see the left-most solid points in the upper panel of Fig. [3]). With data that span 54 years, we have 
shown that at low radio frequencies, the secular decrease is about 0.8% yr~^, and has remained at 
this magnitude for five decades. 

The apparent four modes of oscillation we have found are quite interesting. The DFT of the 
time series displayed in Fig. [J] shows that the oscillations are in fact very consistent with the sum 
of four simple periodic functions. Furthermore, the longer period oscillations seem to be nearly 
1:2, 1:3, and 1:8 resonances of the shortest period oscillation. However, we cannot rule out the 
possibility that the actual variations are somewhat stochastic in nature and that our four-cosine 
model model merely provides an adequate approximation at the time resolution of our data. Still, 
the parameters of the four-cosine fit may hint at possible mechanisms for the generation of the 
oscillations themselves, such as the interaction of shock fronts with both ejecta and the external 
medium; we will leave such issues to be explored by future efforts. 

The ~5-10 yr oscillations also present a problem for using Cas A as a low frequency calibrator. 
While the secular decrease has always made this somewhat problematic, the oscillations worsen 
the problem. Over the 51 years spanned by the data, the four-cosine model fit has an rms of 
about 0.081, corresponding to a variation of 1.24 kJy in the flux density of Cas A from the general 
secular decrease. The model fit also suggests that these variations could be as large as ±2-3 kJy. 
Our four-cosine model does fit the data significantly better than the secular decrease model alone. 
For instance, when the four-cosine model is combined with the secular decrease model, the data is 
reproduced to within ±2.7% on average, with maximum deviations of abut 7%. In contrast, the 
secular decrease model alone does significantly worse, fitting the data within ±5.3% on average with 
deviations as large as 15% or more. However, the likelihood that these four modes of oscillation 
will continue unchanged over a long period of time is unknown. Additionally, over the ~50 year 
period covered by the data and at the data's time resolution, the four-cosine model may simply be 
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a reasonable approximation of a stochastic process that is difficult to predict. Therefore, if there 

are instances where Cas A is needed as a low frequency calibrator, (e.g., to calibrate observations 
of sources which are closer to Cas A on the sky and require phase referencing, or for constraining 
antenna performance for developing northern hemisphere low frequency arrays such as the LWA 
and LOFAR) , one should be aware that the calibration error will be dominated by the limits in our 
ability to predict the flux density of Cas A. 

Overall, it appears that while the secular decrease of Cas A is rather stable at low frequencies 
over 50 years of observations, we have yet to obtain stringent constraints on the characteristics 
and nature of the shorter time-scale fluctuations. While our relatively simple cosine model appears 
to reproduce the data reasonably well, a more complicated set of oscillation modes, or possibly a 
more complex model of stochastic variations may be required to predict the flux density of Cas A 
accurately enough to use it as a low frequency calibrator. Both this and the study of the physics 
behind these fluctuations provide ample justiflcation for a more thorough and consistent effort to 
continually monitor of Cas A at low radio frequencies. 
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of his work with the prototype LWA antennas and for useful feedback. We would also like to thank 
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Table 1. Cas A VLA-based Observations and Flux Densities 













F74 MHz 


Obs. Date(s) 


VLA Proj. ID 


Config. 


(min.) 


(MHz) 


(kJy) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1991-08-25/11-23 


AK272 


A/AB 


131.5 


1.16 


21.1±1.5 


1993-07-08 


AK332 


C 


22.8 


1.16 


22.0±1.0 


1997-05-17 


AR378 


B 


112.3 


1.16 


19.6±0.7 


1998-03-07 


AK461 


A 


21.3 


1.15 


19.3±0.5 


2003-08-02 


AD480 


A 


823.5 


1.15 


15.8±0.8 


2005-02-19 


AP452 


B 


76.5 


1.46 


16.9±0.4 


200G-10-27 


AG 735 


C 


5.3 


1.15 


17.0±2.3 



Note. — Col. (1): Date(s) of the observations. Col. (2): The project ID 
within the VLA archive of the data set used. Col. (3): VLA configuration (s) 
used. Col. (4): The total time spent on source (Cas A). Col. (5): The total 
observing bandwidth used. Col. (6): The integrated 73.8 MHz flux density 
measured using these observations (see §2). 
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Table 2. Four-cosine Fit Parameters 





V 




amp. 






(1) 


(2) 


(3) 


0.074 ± 0.022 


-0.323 ± 0.0025 


-0.10 ±0.085 


0.072 ± 0.025 


-0.111 ±0.0027 


0.15 ±0.097 


0.051 ± 0.024 


-0.044 ± 0.0042 


-0.15 ±0.14 


-0.019 ±0.026 


-0.20 ±0.012 


-0.21 ±0.36 



Note. — Each cosine is of the form Cas A / Cyg 
A = amp cos[27r(z^t ± (/))] where t is time in years; 
the uncertainties for the hsted parameters are the \a 
errors from the covariance matrix computed by the 
fitting routine (see §3.2 and Fig. H]). Note: these 
parameters were derived from a fit to the Cas A to 
Cyg A ratio at 38 MHz. 
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Fig. 1. — For each of the seven epochs for which we have VLA 74 MHz data of Cas A and Cyg A 
(see §2 and Tabled]), the ampUtude of the tapered visibihties for Cas A (see §2) versus u-v radius. 
The data are plotted as black dots with mean values within bins of u-v radius plotted as green dots 
with the standard deviation within each bin represented by the green error bars. Gaussian fits to 
the binned data are plotted in orange. 
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Fig. 2. — The total 74 MHz flux density of Cas A (see §2 for a description of the flux and error 
measurements) as a function of epoch. 
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Fig. 3. — Upper: The ratio of the flux dens ity of Cas A t o tha t of Cyg A at 81.5 MHz for our data 
and other references (see plot legend). The lBaars et al. I (jl977l ) prediction for the secular evolution 



of Cas A at 81.5 MHz, a 1.3% yr~^ decrease, is plotted as a dotted line; a flt to all the data in 
the plot is plotted as a solid line, corresponding to a 0.77% yr~^ decrease. Lower: The r atio of the 
38 MHz flux density of Cas A to that of Cyg A for the data compiled by I Rees I ([l990|) with the 
VLA 74 M Hz data scaled to t he 38 MHz value assuming the spectr al shapes for Cas A and Cyg 
A given by lBaars et al. I (jl977l ). Also i ncluded is the 38 MHz dat a of IVinvaikin I (120061') and scaled 



versions of the 80 MHz BAO data of (jMartirossyan et al. II2002I ). Again, the iBaars et al. I (jl977l ) 
secular decrease predicti on is p l otted as a dotted line. A dashed line is used to represent the 0.8% 
yr~^ decrease derived by lRees I (Il990l ): a fit to all the data is plott e d as a solid line and corresponds 



to a decrease of 0.84% yr . For both the BAO and IVinyaikin I ([2000) data, the data points are 



roughly the sizes of the uncertainties in the Cas A to Cyg A ratio in b oth panels. In both pan els, 
the (scaled) Cas A to Cyga ratio measured with the LWDA antennas bv lHartman et al. I (|2009l . see 
§2) is highlighted with a large open circle. 
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Fig. 4. — Upper: The discrete Fourier transform (DFT) of the data plotted in the lower panel 
of Fig. [3] minus the best-fitting 0.84% yr~^ decrease model (solid black curve) and the DFT of a 
four-cosine model fit to the data evaluated at the same epochs as the data (grey curve). Confidence 
levels computed at each frequency using Monte Carlo simulations are plotted in red, green, and 
blue for the 75%, 90%, and 95% levels (see §3.2). The four largest pairs of peaks which are above 
the 95% confidence level are highlighted in light grey. Lower: The data plotted in the lower panel 
of Fig. [3] minus the 0.84% yr~^ decrease model with a four-cosine model fit plotted as a solid line 
(see Table [2] for the model parameters). The shaded grey area represents the model ztlo", where a is 
the error in the model at a particular epoch based on the covariance matrix of the free parameters 
and the partial derivatives of the model at that epoch. 



